On the Interpretation of "off the edge" Avalanches 
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We establish both experimentally and theoretically the relation between off the edge and internal 
avalanches in a sandpile model, a central issue in the interpretation of most experiments in these 
systems. In BTW simulations and also in the experiments the size distributions of internal avalanches 
show power laws and critical exponents related with the dimension of the system. We show that, 
in a SOC scenario, the distributions of off the edge avalanches do not show power laws but follow 
scaling relations with critical exponents different from their analogous for the internal avalanche 
distributions. 
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Since Bak, Tang and Wiesenfeld (BTW) developed in 
1987 the ideas of self-organized criticality (SOC)P, 01, 
a great amount of research in phenomena as diverse 
as earthquakes, superconducting vortex dynamics, stock 
markets, and ecology d, H, IM 111 has been carried out. 
A sandpile illustrates this concept: the slow addition of 
grains onto a flat surface provokes the growth of a pile 
with slopes around a critical angle adjusted through an 
avalanche mechanism. According to SOC, the avalanches 
should not show any characteristic size or frequency, and 
the distributions of avalanche sizes and durations are ro- 
bust relative to variations of external parameters; i.e., 
the system self-organizes. The result is that the pile will 
show robust power law distributions of avalanche size and 
duration, "1//" power spectra, and finite-size scaling of 
the distribution of internal avalanches, measured as the 
movements of the grains within the totality of the system. 

Although several experiments have been carried out 
trying to find the critical behavior predicted from the 
BTW model in real piles, the results have not shown 
a clear agreement. This is due, in part, to the fact 
that in most experiments only avalanches that involve 
grains abandoning the system (the so-called off the edqe 
avalanches) can be measured. Held and co-workers [3, 
using particles between 1.0—1.3 mm, reported SOC char- 
acteristics when the base of the pile was small enough, 
but when it increased, quasiperiodic large avalanches ap- 
peared. This was later corroborated by Rosendahl et 
aL[3,lal- With a similar setup, but using 3 mm diameter 
spheres, Grumbacher et a/.jl^l and Costello et aLjT^] did 
not find any lost of SOC behavior when the relation base 
diameter /particle size of their piles exceeded the value 
predicted in |3, II, 13 ■ Their piles were also quite sensi- 
tive to the drop height. In "ID" piles, Frette et al.\\3^. 
using rice with grains of different shapes, concluded that 
SOC behavior was attained only for grains with relative 
high aspect ratio, because they were able to decrease in- 
ertial effects, not considered in the theory. In contrast 
with previous experiments, they were able to measure not 
only off the edge avalanches, but those along the surface 
of the pile. Also with elongated rice grains, but in a 3D 
geometry, Aegertcr et al. ^^ found scaling relations mea- 



suring the avalanches as variations in the surface of the 
system. In our previous work on off the edge avalanches 
in "ID" piles |Tj|, the incrtial effects of the 4 mm di- 
ameter spheres used were avoided by introducing strong 
disorder, reached thanks to a base with beads glued with 
random gaps between each other. We claimed that these 
piles show SOC behavior, based on the good critical size 
scaling of the resulting avalanche distributions and on 
the fact that the "active zone" practically involved the 
whole system. As we have seen, none of the experiments 
have measured the avalanches in the whole system, but 
only in a portion of it. In most of them, only the off the 
edge avalanches were considered. 

In the first part of this paper we measure the 
avalanches at several portions of the grid in a two- 
dimensional BTW computer simulation, demonstrating 
that power laws should not be expected in the size dis- 
tributions of off the edge avalanches. After the avalanche 
distributions from the simulations are fitted using an 
appropriate ansatz, we use the resulting equations to 
predict analytically the critical exponents of the corre- 
sponding finite size scalings. The second part reports 
experiments performed in a setup similar to that used 
in [Ijl where the position of the center of all beads are 
obtained after each addition. As a result, both the in- 
ternal and off the edge avalanches are measured, showing 
a non power law for the distributions of the latter, and 
a power law with the same critical exponents expected 
from the 2-D BTW model for the size distributions of 
internal avalanches. 
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FIG. 1: Sites in which avalanches are measured for different 
values of F in a L = 16 grid. 



We perform a BTW computer simulation on a squared 
L X L lattice, following the same rules as in reference [l| . 
The avalanches (defined as the number of sites involved in 
toppling events) were measured in the whole grid (inter- 
nal avalanches) and also in the rows and columns LF — 1 
sites away from the borders, where i^ is a fraction of 
L (see Fig. 1). Then, the avalanche distribution for 
F = 1/L coincides with the off the edge one. The size 
distributions for internal and off the edge avalanches are 
shown in Fig. 2. A cursory inspection of it gives the 
certainty that the distributions of off the edge avalanches 
are not power laws, while the internal avalanches exhibit 
a clear power law behavior over a few decades. However, 
internal avalanches and those corresponding to different 
values of F are just manifestations of the same dynamics, 
so we propose the following ansatz to fit all of their size 
distributions 



p{.s)=pos-'eM-elis/Ny^), 



(1) 



where N is the total number of sites considered for the 
avalanche measurements and e2 = (L/^/2Kf, where 
R{L) is the average distance from the center of the grid 
to all the sites of the zone of measurement \\^. el is 
constant for internal avalanches and el = A F^ /L^ for 
avalanches in the F-dependent sites. The fits to the 2- 
D BTW simulations presented in Figs. 2 and 3 were 
reached with el = 8 ± 1 for the internal avalanches and 
A = 4.0 ± 0.1, 7 = -0.81 ± 0.01 and ^ = 0.28 ± 0.01 for 
those in the F-dependent sites. 

Having reached the equations that fit the tails of all 
the distributions, we can obtain the values of the critical 
exponents for their finite size scalings analytically. If we 
consider a generic slope —a and we write N as kL'^ , where 




FIG. 2: Avalanche size distributions normalized to the total 
number of avalanches for the 2-D BTW model after 5 • 10^ 
additions, for L = 16, 32, 64 and 128. The four curves at the 
right correspond to internal avalanches, while the four curves 
at the left correspond to off the edge avalanches. Black solid 
lines fit Eq. (1) to the data. The scalings based in the ansatz 
p(s,L)L^ = f{sL~'^) are shown in the inset. The critical 
exponents appear in Table 1. 
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FIG. 3: Avalanche size distributions for the F-dependent 
sites (normalized to the total number of avalanches) for F 
equal to 1/2, 6/16, 1/4, 3/16, 1/8, 1/16 and 1/L. Black 
solid lines fit Eq. (1) to the data. 



d is the dimension of the measurement area, we get 



p{s)^Po s-^expi-elis/kLy-") 



(2) 



When the scaling relation p(s, L)L^ = /{sL^") is applied 
to a family of curves (in a log-log plot) they move follow- 
ing a straight hne of slope —a = —fi/v] which suggests 
that if they collapse under scaling relations, it is due to 
the fact that there is an overall power law mechanism 
related with them; therefore, the avalanche distributions 
of F-dependent sites would have the same critical behav- 
ior as the internal ones. In order to analyze that, let us 
define P{S) = p{s)Ll^ and S = sL-", then 



P{S)^Po S-"L''^-°''' exp{-el{SL''-'^/k) 



e2\ 



(3) 



For internal avalanches d = 2, k — 1, and el is constant; 
then we have v = d = 2 and (3 = av . As the slope is 
—a = —1, we get /3 = 2. For the other distributions po ^ 
L^, where 77 = — 0.37±0.01. For constant F distributions 
d = 1 and /c = 4, then 

.(4) 
so i^ = I + ^/e2 and /? — av — rj. We can consider 
a = 1, and calculate the corresponding values for /3; or 
consider (3 — 1.45 (obtained from the values of e2 and 
a for the internal avalanche distributions), and calculate 
the a values. All the critical exponent are shown in Table 
1. For off the edge avalanches d = 1, /c = 4 and F = 1/L, 
then 



PiS) - S'-"L/3-"''+'/ exp(-yl(S'/4)'=2^ 



e2 7--7-?+e2(;/-(i) 



) 

(5) 




FIG. 4: Scalings based on the ansatz p{s,L)L^ = f{sL ") 
for the curves of Fig. 3. The critical exponents are reported 
in Table 1. 

thus, v = l + (f +7)/e2 and /3 = av—rj, with the same as- 
sumption that in the latter case. The inset of Fig. 2 and 
Fig. 4 show all the distributions for the simulations when 
the finite size scalings were applied. The best collapses 
were reached with the critical exponents shown in Table 
1. The match between v values guarantees the quality 
of the fits. A simple inspection of the (3 values from 
the simulation indicates that the correct assumption is 
constant (3. Until here, we have shown that, within the 
BTW model, the distributions of avalanches measured 
in proportional portions of the grid behave similarly as 
in the whole system; i.e., same values of a, and v very 
close to d, with some corrections which depend on how 
far from the center of the system they are. However, 
the distribution of off the edge avalanches do not show 
a power law behavior but collapse when scaling relations 
are applied. The critical exponents do not correspond 
to their analogous from the internal avalanche distribu- 
tions, since their avalanches do not involve sites located 
at constant fraction of the system size, but depend on 
the system size as i^ = 1/L. The relation P/v = 2 is a 
consequence of (s) = 1 for the off the edge avalanches, 
and it is independent from the system dimension. Now 
we present experiments in real piles where these ideas are 
corroborated. 

The experimental setup is similar to our previous one 
[lj| and it consists of an acrylic strip sandwiched between 
two parallel vertical glass plates 5.0 ± 0.2 mm apart from 
each other so that a horizontal surface of 5 x Lmm? 
(where L = 24 cm, 32 cm) was available for the forma- 
tion of a pile of 4.000 ± 0.005 mm diameter steel beads. 
The base consists of a row of the same type of beads glued 
to the surface with random spacings (0, 1, 2 or 3 mm) 
between them. They were dropped one by one from a 
height of 10 ± 2 mm above the apex of the pile. Both lat- 
eral sides of the base were open, so the beads were able to 
fall off the system. The whole setup was computer con- 
trolled in such a way that a digital camera took an image 
of the pile and then a new bead was added only after all 
motion associated to the previous impact had stopped. 
A parallel software found the center of each bead, and 
stored its position. We define an internal avalanche of 
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FIG. 5: Avalanche size distributions (normalized to the total 
number of avalanches) for the experiment. The dashed line 
has a slope equal to -1.15. The scalings based in the ansatz 
p{3,L)L^ = f{sL~'^) are shown in the inset. The critical 
exponents appear in Table 1. 

size s when s beads moved in the whole pile (we define 
"movement" as a displacement of the center of a bead 
not smaller than 1/4 of its diameter when successive im- 
ages are compared) . The experimental equivalent of the 
F sites studied in the simulations, was the measurement 
of avalanches (movements) in different portions of the 
pile (see Fig. 6a). An off the edge avalanche has a size 
s when s beads fell off the pile after a dropping event. 
Each experiment included more than 50 000 avalanches, 
with an average total duration of 250 h. The avalanches 
previous to the attainment of the "steady" average pile 
size were not taken into account for the statistics. 

Fig. 5 displays the size distribution of the internal 
avalanches and of off the edge ones for two system lengths. 
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FIG. 6: (a) Schematic representation of the different zones 
where internal avalanches were measured. The profile is 
the average of two system length profiles (normalized to 
the system length), (b, c, d) Scalings based in the ansatz 
p{s,L)L" = f{sL~^) for the internal avalanche distributions 
for two system length in the different zones. The critical ex- 
ponents appear in Table 1. 



TABLE I: 
Fig. 6a). 



Critical exponents for avalanche size distributions. The experimental analogous of F sites are zones A, B, C (see 



avalanche 


analytic 






simulations 




experiment 




distribution 
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/3(a = l) 


a(/? = 1.45) 


V 


/3 
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internal 


2 


2 


- 


2.00 ±0.02 


2.20 ±0.05 


2.0 ±0.1 


2.3 ±0.1 


off the edge 


0.67 ±0.02 


1.04 ±0.03 


1.61 ±0.07 


0.68 ±0.02 


1.40 ±0.05 


1.4 ±0.1 


2.8 ±0.1 


i^ = l/16 


1.16 ±0.01 


1.53 ±0.02 


0.93 ± 0.03 


1.18 ±0.02 


1.45 ±0.05 






F= 1/8 


1.13 ±0.01 


1.51 ±0.02 


0.96 ± 0.03 


1.11 ±0.02 


1.45 ±0.05 


(C) 2.0 ±0.1 


2.8 ±0.1 


F = 3/16 


1.10 ±0.01 


1.47 ±0.02 


0.98 ± 0.03 


1.07 ±0.02 


1.35 ±0.05 






F= 1/4 


1.08 ±0.01 


1.45 ±0.02 


1.00 ±0.03 


1.06 ±0.02 


1.40 ± 0.05 


(B) 2.0 ±0.1 


2.3 ±0.1 


F = 6/16 


1.05 ±0.01 


1.42 ±0.02 


1.03 ±0.03 


1.04 ±0.02 


1.35 ±0.05 






F= 1/2 


1.04 ±0.01 


1.41 ±0.02 


1.04 ±0.03 


1.04 ±0.02 


1.40 ±0.05 


(A) 2.0 ±0.1 


2.1 ±0.1 



The distributions of internal avalanches show a power law 
regime over two decades with a slope equal to -1.15. The 
critical exponents (also in Table 1) are v = 1 and /? = 2.3, 
very close to those expected from the BTW model in two 
dimensions. This demonstrates that our system is not 
truly one-dimensional, but that the avalanches penetrate 
the bulk of the pile. The distributions of of^ the edges 
avalanches, in agreement with the previous simulations, 
display curves that are not power laws but collapse when 
appropriate scaling relations are applied. In our case it 
happened with v = 1.4 and p = 2.8. The same expo- 
nents were found when the gapran curves from [l4l| were 
re-normalized to the total number of avalanches and the 
scaling relations applied. A value of t^ > 1 indicates that 
the avalanches that leave the pile involve more than one 
dimension, i.e., more than a layer of one bead width on 
the profiles [I^ ■ The exponent values are associated with 
geometrical properties of the system, but they are related 
each other by (3/v = 2. From this view, the power laws 
experimentally found in off the edge avalanches distribu- 
tions in conical piles may be due to their similitude to 
the internal ones as a consequence of the small size of the 
systems. 

The avalanches in different portions of the pile (Fig. 6) 
behave in a similar way to the internal ones for the whole 
system showing the same values of v. The differences in 
area between the zones provoke variations in the values of 
(3. This was corroborated reducing the width of the zone 
B in a 30%. Then both its area and its (i value coincide 
with their analogous for the zone C. All this reaffirms the 
fact that the differences between the critical exponents 
between off the edge and internal avalanches are due to 
the fact that the former ones involve zones that are not 
a fixed proportion of the system length. 

In conclusion, we have explained why several experi- 
ments involving off the edge avalanches do not show the 
expected power laws in their avalanche distributions, but 
they do follow scaling relations. We have also demon- 
strated the existence of SOC behavior in a pile of grains 
provided its components are able to establish a high de- 
gree of inner disorder, needed to spread activity in the 
totality of the system. As a consequence scalings of off 



the edge avalanche size distributions can be safely taken 
as an indication of power law behavior of the internal 
avalanche size distribution in a model sandpile. 
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